Hadron Correlation in Jets on the Near and Away Sides 
of High-pr Triggers in Heavy-ion Collisions 
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Correlation between trigger and associated particles in jets produced on near and away sides of 
high-pT triggers in heavy-ion collisions is studied. Hadronization of jets on both sides is treated 
by thermal-shower and shower-shower recombination. Energy loss of semihard and hard partons 
traversing the nuclear medium is parametrized in a way that renders good fit of the single-particle 
inclusive distributions at all centralities. The associated hadron distribution in the near-side jet can 
be determined showing weak dependence on system size because of trigger bias. The inverse slope 
increases with trigger momentum in agreement with data. The distribution of associated particle 
in the away-side jet is also studied with careful attention given to antitrigger bias that is due to 
the longer path length that the away-side jet recoiling against the trigger jet must propagate in the 
medium to reach the opposite side. Centrality dependence is taken into account after determining a 
realistic probability distribution of the dynamical path length of the parton trajectory within each 
class of centrality. For symmetric dijets with p^ lg = py soc (away) it is shown that the per-trigger 
yield is dominated by tangential jets. For unequal p^ lg , py soc (near) and py soc (away), the yields 
are calculated for various centralities, showing intricate relationship among them. The near-side 
yield agrees with data both in centrality dependence and in py soc (near) distribution. The average 
parton momentum for the recoil jet is shown to be always larger than that of the trigger jet for 
fixed p^ lg and centrality and for any measurable pf? soc (away). With the comprehensive treatment 
of dijet production described here it is possible to answer many questions regarding the behavior of 
partons in the medium under conditions that can be specified on measurable hadron momenta. 

PACS numbers: 25.75.-q, 25.75.Gz, 24.85. +p 
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I. INTRODUCTION 

Recent studies of jet correlation in heavy-ion colli- 
sions at RHIC have generated a wealth of information 
about jet-medium interaction, not only on how the dense 
medium modifies the characteristics of high-p^ jets, but 
also on how intermediate-pr jets affect the medium 
[E S H, Q- Two-particle correlation has been partic- 
ularly effective in revealing the nature of the medium 
response to the passage of a hard or semihard parton 
Q. The discovery of ridge, for example, in the struc- 
ture of the same-side distribution of particles associated 
with a trigger has stimulated intense interest both ex- 
perimentally and theoretically [f|-[l(|. The properties 
of the ridge (in centrality and px dependencies, and in 
baryon/meson ratio) distinguish its origin from that of 
the jet peak that stands above the ridge. Similar distinc- 
tion can be found between the punch-through jet and 
the double-hump peaks on the away side. In this paper 
we study the properties of the associated jets on both 
sides. We calculate not only the px distributions of the 
particles in those jets, but also the fractional energy loss 
of the hard partons traversing the medium toward and 
away from the trigger. As a consequence we can quantify 
the notion of trigger and antitrigger biases. 

It has come to be generally accepted that the 
hadronization process at intermediate px is recombina- 
tion/coalescene I20I l2ll |22| . The approach that we have 
adopted in Ref. 23[ emphasizes the role that shower 



partons play in interpolating the production processes 
from thermal-thermal (TT) recombination at low px to 
shower-shower (SS) recombination at high px, which is 
identical to fragmentation, through the intermediate re- 
gion where thermal-shower (TS) recombination is impor- 
tant. The application of that approach to dihadron corre- 
lation has been considered previously 0, Hf| , but before 
the discovery of ridges. The phenomenology of ridges (or 
ridgeology) has clarified the characteristics of the asso- 
ciated particles on the near side. There are strong in- 
dications that the ridge particles are formed by the re- 
combination of enhanced thermal partons [26j]-[28|]. Thus 
after subtracting out the ridge particles, what remain are 
the jet particles, which being close to the trigger in Acj), 
are due exclusively to TS and SS recombination. With 
the experimental data on the jets being refined, it is now 
appropriate to reexamine the jet correlation problem in 
both the near and away sides. 

Our formalism allows us to study the trigger bias on 
the near side and the antitrigger bias on the away side, 
i.e., higher average initiating parton momentum in or- 
der to allow for more energy loss in traversing longer 
path length to reach the other side. The average trans- 
verse momentum (px) on either side, which is related to 
the inverse slope, turns out to depend sensitively on the 
origin of the partons that hadronize. TS recombination 
has a softer px distribution than SS recombination (or 
fragmentation). The varying mixture of TS and SS com- 
ponents in different px ranges makes the hadronization 
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process an integral part of any procedure to associate 
the characteristics of the px distribution with either the 
trigger or antitrigger effect. We are able to calculate 
the inverse slopes of the associated particles on both the 
near and away sides. Data exist for the former, since the 
ridge contribution has been studied experimentally in de- 
tail and can be subtracted. Our result agrees very well 
with those data. The associated particles in the away- 
side jets are hard to analyze because of the double- hump 
background that is difficult to separate, so our results 
on those jets cannot yet be checked by data. For fixed 
medium suppression in central collisions, we have studied 
various other quantities that are not directly measurable 
in experiments, but they can shed considerable light on 
the trigger and antitrigger effects. 

For realistic nuclear collisions the path lengths of hard 
partons and the quenching effect depend on the point of 
hard scattering in the transverse plane and the azimuthal 
angle of the trajectory. An important part of our study 
is to find a way to describe the variation of both the 
path length and the quenching effect, what may be called 
the dynamical path length, for different collisions within 
each class of centrality. A distribution of that measure 
will play a crucial role in relating theoretical calculations 
to experimental observation at definite intervals of cen- 
trality. Such a distribution has been found in our study 
and is shown to render an excellent reproduction of the 
dependence of the inclusive spectra of pions in the range 
2 < pt < 11 GeV/c at all centralities. With that dis- 
tribution at hand, we can then calculate dihadron cor- 
relation that can realistically describe properties of dijet 
production. We find that the dominance of tangential 
jets emerges naturally as the momenta of trigger and as- 
sociated particle on opposite sides approach each other. 
It is a clear demonstration of the interplay between trig- 
ger and antitrigger effects. 



II. DIHADRON CORRELATION IN THE 
RECOMBINATION MODEL 

We adopt the formalism initiated i n 1 231) for single- 
particle inclusive distribution and in [24| for dihadron 
correlation. Concentrating on only the jet component 
of the associated particles, we ignore TT recombination 
which gives rise to the ridge on the near side and to the 
double-hump on the away side. The medium effect is pa- 
rameterized by an exponential damping factor that de- 
pends on the path length. Our focus is on the momentum 
distribution of the associated particles for trigger-particle 
momentum larger than 4 GeV/c. We restrict our con- 
sideration to midrapidity and study only the transverse 
momentum px- Thus we shall omit the subscript T, and 
use pt to denote trigger momentum, p a for associated 
particle on the near side and p^ for associated particle on 
the away side. Without any subscript, p shall be used 
as a generic symbol for the transverse momentum of any 
hadron. 



In the simplest form the invariant distribution of a pion 

is m 

dNn f dq 1 dq 2 

P—J—= / FqqW,q2)R*{qi,q2,P) (1) 

dp J qi q 2 

where the qq distribution is in general 

F qq { qi ,q 2 )=TT + TS + SS (2) 
and the recombination function (RF) 
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The thermal parton distribution has the form 

T {q x )=q x —^= C Ql e-^ T 
dqi 

so that the thermal pion distribution is exponential 



dN. 



TT 



pdp 
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(3) 



(4) 



(5) 



The quark momentum qi above are just before hadroniza- 
tion at the end of medium expansion. The shower parton 
at that stage is in the vacuum after the hard parton has 
emerged from the medium. Just after hard scattering the 
distribution of the hard parton momentum k of parton 
type i, while still in the medium, is given by 



dN, 



hard 



kdkdy 



fi(k), 



(6) 
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whose specific properties are given in the next section. 
After propagating through the medium, the parton loses 
momentum in a way that we describe by the function 
G(q, k, t) where t denotes the distance the parton travels 
to reach the surface, and q is the momentum of that par- 
ton at the surface. We discuss G(q, k,t) below presently. 
The parton distribution in q after averaging over all k 
and t is 



Fi{q) 



L dt r°° 

- dkkfi(k)G(q,k,t), 

L Jkn 



(7) 



where in calculation we set the lower limit fco at 3 GeV, 
below which the parton distribution fi(k) is not known 
reliably. L is the average maximum length of the system 
that the hard parton traverses. In the limit L — * 0, we 
should recover the parton distribution for pp collision. 

For the degradation factor G(q, k, t) due to energy loss, 
there is a rich literature on the subject studied by various 
methods. Two articles reviewing the subject are Refs. 



2S,|30(. The quenching factor Q(jp) determined in Ref. 

| increases with p, a property at very high energy not 
found at RHIC. In the opacity expansion approach [32| 
the energy loss is found to depend on the path length as 
AE oc L 2 ~ a , where a = 1 for one-dimensional expansion. 
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The effective quark energy ioss with detailed balance be- 
tween induced gluon emission and absorption taken into 
account has the form for a 1-d expanding medium [33[ 

=e (E/»-1.6) 1 */(7.5 + E/ri, (8) 



dL/ la 

which is essentially oc E for 4 < E < 12 GeV. A reason- 
able summary of these properties is 



E 1 ' 



(9) 



whose implication for the relationship between q and k 
in Eq. is that 



k — q = k(3t. 



(10) 



For t not infinitesimal, we exponentiate Eq. (|10p and get 

q = ke- pt . (11) 

Fluctuation from this relationship is undoubtedly possi- 
ble, but we shall take the simple form 



G(q,k,t) =q6{q-ke- l3t ) 



(12) 



as an adequate approximation of the complicated pro- 
cesses involved in the parton-medium interaction. The 
justification for Eq. (fl"2"|) is to be found in the degree 
to which the inclusive cross section can be reproduced 
in our description of hadron production at intermediate 
and high pt- 

Using Eq. (O, the integration over t in Eq. ([7]) can 
readily be carried out, giving 



1 f qefi 



dkkfi{k) 



(13) 



for q > fco- The lower limit of the above integration cor- 
responds to t = 0, i.e., when the hard scattering occurs 
at the surface, while the upper limit corresponds to the 
hard-scattering point being on the far side so that k is a 
factor e^ L larger than q. Equation (| 13[) exhibits the nu- 
clear effect in changing fi(k) to Fi(q) with j3L being the 
explicit medium factor, while fi{k) contains the hidden 
modification of the parton distributions in the nucleus 
due to shadowing, etc. [HI Clearly, as f3L — > 0, Fi(q) 
becomes directly related to fi(q) appropriately extrapo- 
lated to pp collision. 

Using Si to denote the matrix of shower parton distri- 
butions (SPDs ) th at are calculable from the fragmenta- 
tion functions [351 ], we can determine the distribution of 
shower partons in a heavy-ion collision by 



S(qi) = J ^-FMSjiqi/q). 



(14) 



The TS contribution to the inclusive pion distribution is 
then, following Eqs. Q, ©, © and 



dN? s 
pdp 



v J q 



(15) 



with the RF absorbed in the compound notation for the 
TS term in the integrand: 

TS(?,P) = / ^S> (^j J dq 2 C ] e~^/ T R fj (q 1 ,q 2 ,pU6) 

where for every hard-scattered parton of type i the 
shower parton of type j is paired with a thermal par- 
ton of type j for recombination in forming a pion. For 
the SS component we can use the fragmentation function 
D(z) and write 



dN. 



ss 



1 



pdp 



P 
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(17) 



The overall pion inclusive distribution, including the TT 
contribution as given in Eq. (0, is thus 



pdp 



C 2 



-p/T 
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TS(q,p) + ^DJ 

q 



dq 

p 
<i 



Fi(q) 



(18) 



For the dihadron correlation on the same side we con- 
sider the trigger-momentum p t to be greater than 4 
GeV/c and calculate the associated particle distribution 
in the approximation that the TT contribution to the 
trigger and jet can be neglected. We then have for p a on 
the near side associated with p t 



1 



PtPadptdpa (PtPa, 



2 Z-^i 



[ "Mia) 

J q 



TS(q lPt ) + ^D? 

q 



TS(q-pt,p a ) 



-fS(g-p a ,p t )-AM- 

q \q 



PtPa m 




(19) 



where the dihadron fragmentation function 1)2(^1,^2) is 
assumed to have the symmetrized form 



D 2 (z ll z 2 ) 



-D 



D( Zl )D 

Z\ 



Z-2 



l-z 2 



l-zi 
D(z 2 ) 



(20) 



The near-side yield per trigger for trigger momentum in 
a narrow range Ap t around pt is 



(j>t,Pa) 



1 rfiNU 

Ntrlg Padpa ' 

dpt , 

Ap t Padptdpa 



, dN v 
dpt-; — 

A Pt dPt 



(21) 



where dN^/dpt is the trigger pion distribution that ex- 
cludes the TT component of the inclusive distribution 
given in Eq. ([T5jl. 
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For an associated particle on the away side relative to 
the trigger we must consider the recoil hard parton that 
propagates a distance L — t in the opposite direction, so 
Eq. ([7]) should be revised to contain another parton with 
momentum q' exiting on the away side, having the form 



% I dkkMk)G(q,k,t)G(q',k,L-t) 



L Jk 



;5L 



= Jlj dkkfr(k)qq'5(qq' - k 2 e^ L ), (22) 

where the recoil parton has momentum k under the as- 
sumption of negligible initial fcy of the beam partons. 
However, for reasons that will become clear later we label 
it by k' to be distinguished from k of the trigger parton, 
when clarity is needed. To help with the visualization of 
the various momentum variables in the problem, we give 
a sketch in Fig. 1 of the parton (red) and hadron (blue) 
momentum vectors with the trigger being on the right 
side. With due caution in the interpretation of momen- 
tum vectors drawn in coordinate space, momenta on the 
near side (q,pt and p a ) are depicted to originate from the 
surface on the right, while the momenta on the away side 
(q' and pb) are pointed from the surface on the left. 




q « 



FIG. 1: (Color online) A sketch of momentum vectors of par- 
tons (in red) and hadrons (in blue) with near side being on 
the right and away side on the left. 

The recoil parton of type i' need not be linked to the 
parton type i of the trigger jet, since beam partons can be 
of any type. However, to limit the problem to a manage- 
able size we make the simplifying assumption that %' = i 
if they are quarks and i' = i if they are gluons (thus tan- 
tamount to assuming dominance by g+g scattering), and 
we shall only calculate identical pions on the two sides. 
The dipion distribution for particle with pb on the away 
side is thus 



1 



j J q q 



ptpbdptdpb 



iptPb) 



q \q 

TS (q ,p b ) + —DJ, — 

q \q 



Carrying out the integration over k yields 



(23) 



d/NL 



PtPbdptdpb 2l3Lpfpl 



E / dq r dq'fi U/tfrfZ 



TS((?,_pt) 
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Pb m 
7° v 


( Pb 





(24) 



where q' = Ma,x(qe~ f3L ,pb). 

We note that the limits of integration of q' in Eq. 
reveals the medium effect in the following sense. If the 
hard scattering occurs at the near-side surface, then the 
recoil parton (having k' = k = q) must travel a dis- 
tance L to emerge on the away side with momentum 
q' = q'o = qe~^ L . On the other hand, if the hard scat- 
tering occurs at the away-side surface, then q must be 
ke~^ L , so q' — k = qe l3L . Thus the integration over q' 
reflects the integration over all points t in the medium 
where the hard parton is created. 

The away-side yield per trigger is, analogous to Eq. 
(EH), 



Y^{p u pb) 



1 dN„ 



N, 



-(j>t,Pb) 



trig Pbdpb 

dN 
dpt , , 
Apt Pbdptdpb 



d Pt -^(25) 
A Pt dp t 



In our calculation we can, of course, take the theoretical 
limit Ap t — * 0. 



III. MODEL INPUTS 

We list here all the inputs to the model that we use to 
perform our calculation. They are all taken from previous 
work without any parameters to adjust, except for (3L 
that is introduced here in lieu of an average suppression 
factor used earlier. 

The hard-scattered parton distribution fi(k) is taken 
from Ref. [3~4| . which uses the parametrization 



fi(k) = K 



(l + k/B) a 



(26) 



with K = 2.5 and the parameters A, B, and a tabulated 
for each parton type i and for Au+Au collisions at RHIC 
with shadowing taken into account. The sum J^. will be 
performed over i = g,u, d, s, u, d, s. For the thermal par- 
tons the values of C and T in Eq. ([5]) for 0-10% centrality 
are [H 



C = 23.2 GeW 



T = 0.317 GeV. 



(27) 



Their centrality dependence are given in Ref. [25j. The 
shower parton distributions are described in Refs. [H, 
35]. For z in S\{z) very small, the distributions are not 
reliable, so we cut off the low-pT contribution to the TS 
component by a factor 1 — exp(— O.bpx), which has no 
effect on our result for intermediate and high px- For the 
fragmentation function D(z) we use the parametrization 
in Ref. [36| . from which the shower parton distributions 



5 



were derived [351 ] . Since we now consider higher px than 
before, the Q 2 dependence of D(z, Q 2 ) will be included 
by setting Q 2 = p\. 

With these inputs specified there are no more free pa- 
rameters to adjust, except the suppression factor quanti- 
fied by PL, which will be determined below by fitting the 
overall single-pion inclusive distribution. The properties 
of the dihadron correlations on both the near and away 
sides can then be calculated without unknown parame- 
ters. 



IV. NEAR-SIDE CORRELATION AND 
TRIGGER BIAS 

We first calculate the pion inclusive distribution using 
Eq. (fl~8|) and compare the result to the data in Fig. 2. 
What we have calculated is dN^ / prdpr at midrapidity 
averaged over all </>, while the data are for dN n o /2irpTdpT 
integrated over all <f> [13], both for 0-10% centrality in 
Au+Au collision at y/s — 200 GeV. The value 



PL = 2.9 



(28) 



has been used to fit the data for 2 < pr < 13 GeV/c. 
Since the suppression factor involving PL enters Eq. Q18p 
only through Fi(q) given in Eq. (fT3"]) , the excellent fit 
in Fig. 2 over such a wide range of pt requires a high 
degree of coordination among the three components of 
recombination, and therefore is not a trivial result from 
varying one quantity, PL. The agreement with data is a 
confirmation of the soundness of the model for the range 
of pt considered. 




4 6 ,J3 10 12 14 
p T (GeV/c) 



FIG. 2: Inclusive distribution of 7r° as calculated from Eq. 
(fig)) . The data are from Ref. [37[ . 



Using PL given in Eq. (f2"5| . we can now calculate the 
dihadron correlation on the near side for trigger momen- 
tum at pt > 4 GeV/c and an associated particle in the jet 
with momentum in the range 2 < p a < 4 GeV/c. We use 
Eq. ([TO)) to study the irir correlation without TT recom- 
bination, which we have assigned to the ridge. The result 
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FIG. 3: (Color online) Distribution of associated pion (p a ) in 
near-side jet for six values of pion trigger momentum (pt) in 
GeV/c. 



for the yield per trigger is shown in Fig. 3 for 6 values of 
Pt- It is evident that the p a spectrum becomes slightly 
harder, as pt increases. The effective inverse slope T a de- 
termined in the range 2 < p a < 4 GeV/c is shown in Fig. 
4. The data in that figure are from Refs. @, 0] for all 
charged hadrons. Although our result is for pions only, 
we expect that the contributions from the other charged 
hadrons are not as important in the jet as they are in the 
ridge. Thus the general agreement of our result with the 
data may be regarded as supportive of our description of 
the physics that generates the dihadron correlation. 




5 6 7 8 
p t (GeV/c) 



FIG. 4: Inverse slope of associated pion in near-side jet. Data 
are from Refs. @,0] determined from 2 < p^ BOC < 4 GeV/c. 



The results obtained so far average over all possible 
points of hard scattering, as indicated by the integration 
over t in Eq. ([7]). However, we know that with pL = 2.9 
the contributions from the points on the far side of the 
medium are more suppressed compared to those from the 
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nearer points. That is the trigger bias in heavy-ion col- 
lision. We can quantify that effect by calculating the 
average of exp(— flt). Such an averaging process is feasi- 
ble by using Eqs. (J7J) and (TT5)) . where we insert exp(— f3t) 
in Eq. ([7]) before integration over t. It is not necessary to 
know separately from (3L because of the structure of 
G(q,k,t) that demands e~P* = q/k. With the near-side 
suppression factor defined as 

rncarOPT) = (e^*), (29) 

which is also (q/k), we show our calculated result in Fig. 
5. Evidently, it saturates at 0.85. Thus on average only 
15% of the parton energy is lost to the medium when 
Pt is high, but more at lower px- This result is roughly 
independent of the medium size L, and is a feature of the 
trigger bias. That is, if the point of origin is allowed to 
vary, the detected hadrons are dominantly due to partons 
created near the surface and losing only a small fraction 
of the energy. 

1i ■ ■ ■ ■ ■ 1 



0.9 




c 

0.7 

062 4 6 p t( « 12 14 

FIG. 5: Near-side suppression factor for which the value of 
1 means no suppression. The averaging is done over single- 
paticle distribution of pion momentum pr- 




FIG. 6: Average values of parton momenta k (at point of hard 
scattering) and q (at the near-side surface) for pions detected 
at pt- 
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FIG. 7: Average value of fit as a function of pion momentum 
Pt- 



In Fig. 6 we show both (k) (in solid line) and (q) (in 
dashed line) as functions of pr- Their ratio (q) / (k) is 
not exactly (q/k). The two lines provide insight into 
how hard the hard scattering has to be in order to give 
rise to a pion at px- Note that (k) is approximately 
2pT throughout the range, while (q) « 1.6 px for pr > 3 
GeV/c where both TS and SS components of recombi- 
nation are important. Since (exp(— (3t)) ^ exp(— /3<), we 
have calculated (flt) shown in Fig. 7. It is an estimate of 
AE/E according to Eq. ©, although TS recombination 
renders the connection with parton energy loss impre- 
cise. For pt > 4 GeV/c the average ((3t) is between 0.18 
and 0.2. That is to be compared to (3L — 2.9, implying 
(t)/L » 0.065. Thus the result suggests that the thick- 
ness of the layer near the surface where hard partons 
are created is roughly 13% of L. That is a quantitative 
statement about the trigger bias to the extent that we 
can calculate without taking into account such details as 
nuclear geometry, tangential jets, etc. 



V. AWAY-SIDE CORRELATION AND 
ANTITRIGGER BIAS 

As we consider the correlation between jets on oppo- 
site sides, we first fix f3L = 2.9, which corresponds to 
a slab of nuclear medium with fixed thickness. That is 
not a realistic nuclear medium, whose thickness depends 
on the transverse distance from the center. The value of 
f3L determined in Eq. (128|) corresponds to the effective 
thickness in fitting the single-particle distribution shown 
in Fig. 2. Due to trigger bias the correlation between 
particles on the same side is mostly independent of that 
thickness, as we have seen in the preceding section. Now, 
as we go to dijet correlation on the two sides, it makes 
a big difference whether the medium has varying thick- 
ness. In order to illuminate the nature of the away-side 
correlation and antitrigger bias, we first consider in this 
section the simplest scenario of fixed thickness. After be- 
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coming familiar with the issues involved, we then extend 
our study to the case of realistic nuclear medium in the 
next section. 




P b (GeV/c) 

FIG. 8: (Color online) Distribution of associated pion (pf,) in 
the away-side jet for six values of pion trigger momentum (pt) 
in GeV/c. 
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FIG. 9: Inverse slope of associated pion in away-side jet. 

For the away-side yield per trigger we use Eqs. (f24|) 
and p5|) to calculate Y£™ ay (pt,Pb) as functions of pb for 
six values of pt- The results are shown in Fig. 8. As with 
near-side yield, the away-side per-trigger yield increases 
with pt, not simply because the corresponding hard par- 
ton k is forced to be higher, but also because the number 
of triggers is lower. However, the spectrum does not be- 
comes harder at larger p t , for a reason to be discussed 
later. Figure 9 shows the inverse slope T\, determined 
in the range 2 < pb < 4 GeV/c, exhibiting only a mild 
decrease of ~ 10% over the range of p t - Since it is a 
property of the jet yield, there are no suitable data to 
compare with our result. PHENIX has extensive data 
on dihadron correlation however, on the away side 
the division between head and shoulder regions is done 



in terms of cuts in A<f>, with the consequence that a di- 
rect relationship between the yield in the head region and 
the jet yield calculated here cannot easily be established. 
Data on inclusive 7 have been analyzed for correlated 
hadrons, using a Gaussian description for punch-through 
jets on the away side [38j ] . but no px distribution has 
been shown. 

To learn about the medium effect on the away-side jet, 
we study the suppression factor defined as 

r awa y(p t! P6) = (exp[-/3(L - *)]). (30) 

It should be recognized that whereas r noar (pT) in- 
volves an average over the single-particle distribution, 
r away (p t ,pb) requires for the averaging process the 
opposite-side dihadron correlation that depends on p t 
and pb- Our results are shown in Fig. 10 as functions 
of pt for four values of pt- As expected, the suppres- 
sion is far more severe on the away side than on the near 
side. The higher the trigger momentum p t , the closer is 
the hard-scattering point to the surface on the near side 
due to trigger bias, and we see that the more severe is 
the suppression on the away side. That is a property of 
antitrigger bias because of the longer path length that 
the recoil parton must travel in the medium. Whereas 
r n ear(pr) increases with px, here r away (p t , pb) decreases 
with pt, but increases with pb for a fixed pt- That is be- 
cause higher pb requires higher hard-parton k or shorter 
L — t. Fixing p t does not fix k, as the value of (k) in Fig. 
6 is determined without any extra condition. Now, with 
higher pb, k must be higher, as well as t must be larger, re- 
sulting in increasing (exp[— j3{L— t)]). In other words, the 
condition of higher pb favors the hard-scattering points 
closer to the away-side jet in order to reduce suppression. 



p =4 GeV/c 
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FIG. 10: Away-side suppression factor determined by aver- 
aging over opposite-side tttv correlation function with trigger 
momentum at p t and associated-particle momentum at pb- 

In a hard scattering process the two outgoing hard par- 
tons have equal and opposite momenta, if we ignore the 
transverse momenta hr of the initial partons. However, 
the averages of the two hard parton momenta may dif- 
fer, depending on what observables are held fixed. That 
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FIG. 11: Average value of recoil parton momentum k for 
various values of trigger p t and associated-particle pb in the 
away- side jet. 



is, if k! is the recoil momentum, opposite to k that gen- 
erates the trigger, we have k' = k event-by-event, but 
(k')(pt,Pb) may well be different from (k)(pT)- We have 
seen how (k)(pT) depends on px in Fig. 6. We now show 
in Fig. 11 the dependence of (k')(pt,Pb) on pb for four 
values oi pt- Evidently, (k'} is much larger than (k) for 
all values of pb and pt- This is the essence of antitrig- 
ger bias. The condition of having a hadron on the away 
side among the triggered events gives higher weight to 
the larger k' processes in the averaging. That is why (k') 
increases with both p t and pb ■ 

As pb increases, (k') must increase in order to provide 
enough (q 1 ) that can accommodate the larger pb- Figure 
12 shows (q'){pt,Pb), which exhibits a mild dependence 
on p t , but rises almost linearly with pb above 3 GeV/c 
as ~ (1.5 -j- 2)pb- That dependence on pb is roughly the 
same as the dependence of (q) on px shown in Fig. 6, 
as it should, since the hadronization processes are simi- 
lar. The weak dependence on p t implies that as pt is in- 
creased, the point of hard scattering is pulled to the near 
side at the same time as the scattered-parton momentum 
k is increased (see Fig. 6), so their opposite effects on the 
away side q' are nearly canceled. With (q r ) not changing 
much in the range of pt probed, the jet yield on the away 
side remains roughly the same, but the per-trigger yield 
increases due to the decrease of the number of triggers 
at higher p t ; that property of Y£™ ay (p t ,pb) is shown in 
Fig. 8. Since (q 1 ) is insensitive to pt, the shape of the 
Pb distribution should therefore also be insensitive to pt , 
and that is confirmed by the property of X& in Fig. 9. 

The values of (q 1 ) are smaller than (k') because of the 
longer path length for the recoil parton to reach the away 
side. We have already a hint of that in Fig. 10, since 
r away (Pt,Pb) is also (q' /k')(p t ,Pb), owing to the 5 function 
in G(q', k',L-t), despite the fact that (q'/k') (q')/(k') 
exactly. Nevertheless, at p t — Pb — 4 GeV/c, both 
(q')/(k') and {q'/k 1 ) are approximately 0.24. 




FIG. 12: Average value of parton momentum at the away- 
side surface for trigger momentum pt and associated particle 
momentum pt- 

Let us now make a more careful comparison between 
(k) and (k'), bearing in mind the difference in the prob- 
ability distributions used in the averaging of k and kl . If 
we set pt in Fig. 6 to p t in Fig. 11 and consider pb — pt in 
the latter, then we find (k 1 ) = 4.7 (k) forp-r — p t — pb = 4 
GeV/c, (dropping to 4.2(fc) at p t = 6 GeV/c). Thus 
the condition of having an associated particle on the 
away side with equal momentum as the trigger elimi- 
nates the trigger bias and raises (k 1 ) to approximately 
9 pt- The implication is that the average location of the 
hard-scattering point is in the middle of the medium (due 
to the symmetry of the two sides) and a large fraction of 
the parton energy is lost before exiting on either side. 
That fraction is 76% as we have learned from (q') j (k 1 ) 
above. 

There are properties of the suppression factor that are 
noteworthy at the symmetry point p t — Pb- We see in 
Fig. 10 that F away (p t ,pb) appears to be constant when 
Pt = Pb is changed from 4 to 6 GeV/c. Let us then define 

r 0) = Laway(Pt,Pb), P = Pt = Pb, (31) 

and calculate its p dependence. The result is shown in 
Fig. 13. It is indeed constant with the value T(p) = 0.24. 
As we have stated above, it follows from the 6 function in 
G(q',k',L-t) that r away (p t ,p b ) = (q'/k')(p t , Pb ). Thus 
for p t — Pb, we have 

T(p) = (q'/k') = (q/k) = e-^ L ' 2 , (32) 

the last equality being the consequence of identifying Eqs. 
(|29|) and (|30| at the symmetry point. Putting (3L = 
2.9 in Eq. (J32J), one finds that T(p) = 0.235, which is 
essentially the value 0.24 determined in Fig. 13. Thus we 
have consistency. 

It is also of interest to compare (q) and (q 1 ) for cor- 
related particles on both sides. The result on (q') / (q) is 
shown in Fig. 14. Note that the ratio is 1 at the sym- 
metry points pt = Pb at 4 and 6 GeV/c. For fixed pt, 
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FIG. 13: Away-side suppression factor for symmetric mo- 
menta on the two sides: p — p t — Pb- 
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FIG. 14: The ratio (q')/{q) as a function of pt and pb- It has 
the value 1 at p t = Pb, shown explicitly at 4 and 6 GeV/c. 

the ratio increases with pb, and, of course, for fixed pb it 
decreases with p t , since the inverse ratio, (q) / (q') , should 
increase. This feature essentially describes how the hard- 
scattering point moves from one side of the midpoint to 
the other side, as pb is changed from below p t to above 
p t . That is clearly a consequence of the counteracting 
effects of trigger and antitrigger biases. 

VI. CENTRALITY DEPENDENCE 

In the preceding section we have studied the proper- 
ties of dihadron correlation on opposite sides for a nuclear 
slab of fixed /3L = 2.9, which is the average value of the 
quenching parameter in central collision. Before consid- 
ering other centralities we must first establish a scheme 
to treat the realistic nuclear medium that has varying 
transverse width even for a fixed centrality, depending 
on the section of the overlap that the hard parton tra- 
verses. As a start we can formally treat PL as a variable 
and study how the yield depends on it. 
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FIG. 15: Yield per trigger at the symmetry point p = p t — Pb 
for (a) fixed £, and (b) fixed p. 

Let us define the per-trigger yield at the symmetry 
point 

Y(p,Z)=Y£*(p t ,p b ,!3L), (33) 

where p — p t — pb and £ = [3L now treated as a variable. 
Using Eqs. (f24l) and (|25|) . we can calculate Y(p, £) with 
the results shown in Fig. 15(a) and (b) for fixed £ and 
p, respectively. The decrease with p for fixed f is more 
gentle than in Fig. 8 for unsymmctrical pt and pb- The 
decrease with £ for fixed p is exponential for £ > 0.5, 
approximately as e -3 ' 8 ^, which is roughly what one ex- 
pects from the j3L dependence that one sees in Eq. (|24| , 
remembering that f%(k) behaves in a power law as indi- 
cated in Eq. (|26[) with a ~ 7.7 -j- 8.7. What we gain from 
Fig. 15 is a quantification of the picture that we already 
have in the increase of yield when the thickness of the 
nuclear medium is decreased. A corollary to that picture 
is that the hard parton momentum k' need not be much 
larger than p, when £ is smaller. That is shown in Fig. 
16, where (k') decreases by nearly a factor of 3 when £ is 
decreased from 3 to 0.5. 

The above consideration is merely a preview of what 
one should anticipate when we treat the medium realis- 
tically and change the centrality. What we need first is 
a distribution of £ for a fixed centrality. Since quenching 
characteristic is involved in PL, it is not just a geometri- 
cal problem of determining the path length in the elliptic 
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FIG. 16: Average value of recoil parton momentum as a func- 
tion of p and £. 
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overlap. The quenching effect depends on the local nu- 
clear density and the location and orientation of the par- 
ton trajectory, so £ is a measure of the dynamical path 
length. Obviously, it is a very complicated problem for 
which no reliable solution is known. We shall approach it 
by first determining the average (3L for single-particle in- 
clusive distribution at each centrality c, as we have done 
in Fig. 2 for 0-10% centrality. We use c to denote the % 
centrality so that c = 0.1 means 10% centrality, for exam- 
ple. We then construct a probability distribution P(£, c) 
such that the average £(c) can fit the average (3L as a 
function of c. With P(£, c) at hand, it is then possible to 
calculate the yield per trigger for any centrality. 

We start by revisiting Sec. 4 and the beginning of Sec. 
5, but now consider all centralities c. For thermal par- 
tons the values of C and T in Eq. (4) as functions of c 
are given in Ref. [25| , For hard partons, fi(k) is scaled 
by N co i\. We use Eq. (fig)) to calculate the single-pion pr 
distribution and obtain a good fit of the data, shown in 
Fig. 17, by adjusting (3L(c) at each c. The data are from 
PHENIX [13 for c = 0.05, 0.15, 0.86 at intervals of 
0.1. The fits are remarkably good for all centralities. The 
resulting (3L(c) are shown by the solid dots in Fig. 18, 
which will serve as the key link between the centrality de- 
pendence of realistic nuclear collisions and the modeling 
of the quenching probability P(£, c) at each c. 

With the aim of fitting (3L(c) in mind, it is sufficient 
to use a simple form for P(£, c) that contains the basic 
features of noncentral collisions, namely: at any fixed c, 
P(£, c) should have a maximum between the two ends 
of £, with the location of the maximum decreasing with 
increasing c. We adopt the form 



(34) 



where N normalizes the total probabilty to 1, and £o arL d 
a are two parameters. We find that with 



£o = 5.42, 



a = 15.2 



(35) 



we get the average £(c) that fits f3L(c) very well, as shown 
by the solid line in Fig. 18. The distribution P(£, c) it- 



FIG. 17: Inclusive distribution of tv° for all centralities rang- 
ing from 0-10% (top) to 80-92% (bottom) in 10% steps, each 
displaced by a factor of 0.2. The data are from Ref. [37|. The 
curves are calculated from using Eq. (|18[l with (3L{c) adjusted 
to fit and shown as dots in Fig. 18. 




FIG. 18: The dots are the values of /3L(c) used to fit the 
inclusive distributions in Fig. 17. The solid line is the average 
^(c) from the £ distribution P(£, c) given in Eq. (|34[) . 




FIG. 19: The distribution of the dynamical path length £ for 
six values of centrality c in steps of 0.1. 
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self is shown in Fig. 19 that exhibits the decrease of the 
maxium with increasing c. In view of the difficulty of 
deriving j3L{c) from first principles, let alone P(£,c), we 
regard Eqs. (|3~i|) and ([55)) as being totally satisfactory 
for the description of how the path-dependent quenching 
parameter varies among the collisions within each class 
of centrality. 

With P(£, c) thus obtained, we can now return to the 
per-trigger yield Y(p, £) at the symmetry point defined 
in Eq. ([55]) . To determine the yield at a fixed centrality 
it is not simply a matter of averaging Y(p, £) over all £, 
using P(£, c) as the weighting factor at each £. Y(p, £) 
that is shown in Fig. 15 is obtained for centrality being 
held at c = 0.05, while £ is varied. We must redo the 
calculation for Y^(j> t ,pb) using Eq. (24|) and ([25]), but 
now include also the dependencies of C, T and fi(k) on 
c, discussed above. That is, we define 



Ptdpt 



^(£,c) 



PtPbdptdpb' 
dN^(c) /dN w (c) 



ptp b dp t dp b 
Y^[p t ,p b ,c) 

and calculate 

Y(p,c)=Y™*y(p= Pt =p b ,c 



PtPbdptdpb I ptdpt 



(36) 
(37) 
(38) 

(39) 



The results are shown in Fig. 20 for four centralities. 
Note that the dependence on c is not as drastic as the 
dependence of Y(p,£) on £ in Fig. 15(a), which shows 
the dominance of £ = 1 over £ = 3, so upon averaging 
over £ at each c the small £ contribution is always more 
important at any c. The per-trigger yield rises with c be- 
cause of reduced medium suppression. The p dependence 
appears to be universal. 
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FIG. 20: Yield per trigger at the symmetric point p = p t — Pb 
for four values of centrality. 

Related to the mild dependence on c in Fig. 20, we 
can investigate the properties of (k'). Recall from Fig. 
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FIG. 21: Average value of the initial hard parton momentum 
directed at the away side for hadron momenta p — p t — Pb for 
four values of centrality. 
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FIG. 22: Average value of the parton momentum at the away- 
side surface for hadron momenta p — p t — Pb for four values 
of centrality. 



11 that for fixed /3L = 2.9, (k') is much larger than p t or 
Pb, the phenomenon referred to as a feature of antitrig- 
ger bias. We have also seen that at the symmetry point 
p — Pt — Pb the values of {k')(p, £) decrease significantly 
at lower £, shown in Fig. 16. To calculate (k'}(p,c) for 
different c we again cannot simply average (k')(p, £) over 
£ using P(£, c) as weight, since the normalization factor 
of (k 1 ) (p, £) must also be averaged over £ separately. The 
result for (k')(p,c), shown in Fig. 21, exhibits essentially 
no dependence on c. The magnitude is approximately 
2.5p, which is much lower than £) in Fig. 16 and 

more like the near-side (k)(pr) in Fig. 6. Furthermore, 
(q')(p, c) can also be calculated in the same manner with 
similar result shown in Fig. 22. The small difference is 
that whereas (k'){p, c) decreases slightly with c, (q')(p, c) 
increases imperceptibly. The ratio (q') /(k') is seen in Fig. 
23 to be nearly constant at around 0.8, increasing about 
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4% when c changes from 0.05 to 0.35. The value of that 
ratio is roughly the same as the value of r near (p^) at com- 
parable pt in Fig. 5, which corresponds to (q/k) on the 
near side without the requirement of a recoil jet. Thus 
the medium degrades the parton momentum from k' to 
q' on the away side by about the same degree as from k 
to q on the near side, and the degree of suppression is es- 
sentially independent of centrality. The inescapable con- 
clusion is then that when symmetric back-to-back hadron 
momenta (j> t = pt,) are required, the dijets that give rise 
to them are due to hard partons created very near the 
surface on both sides so that they suffer minimal energy 
loss as they propagate in opposite directions through the 
rim of the nuclear medium. That means they must be 
tangential jets. This is a remarkable result that emerges 
from the calculation, and is consistent with the dijet+1 
correlation data [39j in which no ridge is found and whose 

2/3 

N p ' aIt dependence suggests that they are generated near 
the surface, i.e., tangential jets. 



for 2 < p a < 4 GeV/c. The data in Fig. 24 are for 
3 < Pr ig < 4 GeV/c and pf? soc > 2 GeV/c [40]. The 
agreement between our result and the data is remark- 
ably good both in magnitude and in c dependence. The 
magnitude of the integrated yield is sensitively dependent 
on the lower limit of integration in p ai so other data on 
centrality dep endence with different lower limits, such 
as in [18|, |4l|, cannot be compared with the black line 
in Fig. 24, although the rough insensitivity to N paxt is 
seen irrespective of the cut in p|f soc . The approximate 
independence on centrality is a manifestation of the trig- 
ger bias, as we have already noted in Sec. 4 that the 
hard-scattering point is in a layer roughly 13% of L in- 
side the near-side surface and is insensitive to how large 
the main body of the medium is. However, the thermal 
and shower partons have different dependencies on c and 
the decrease of TS recombination with increasing c can- 
cels the increase of SS recombination with c so that their 
sum results in approximate independence on c. 
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FIG. 23: The ratio {q')/{k') at the symmetry point p = p t 
Pt for four values of centrality. 



VII. YIELDS AT UNEQUAL TRIGGER AND 
ASSOCIATED PARTICLE MOMENTA 

Having studied in the previous section on how the yield 
at the symmetry point p = pt = Pb behaves at differ- 
ent centralities, we can finally investigate the properties 
at asymmetric points and appreciate the significance of 
small variations. We first consider the per-trigger yield 
of an associated particle on the near side at centrality 
c. The basic formula is as in Eq. (f2Tj) . except that both 
the numerator and denominator must be averaged over 
P(£,c) separately, similar to Eqs. (|3"o]) - (f55|) . The cen- 
trality dependence of the result is shown in Fig. 24 for 
Pt = 4 and 6 GeV/c and p a = 2 and 4 GeV/c. The near- 
side yield is nearly constant in c, and decreases with p a 
for a fixed pt, but increases with p t for a fixed p a - The 
solid lines in that figure represent the integrated results 
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FIG. 24: (Color online) Yield per trigger in the near-side jet 
as functions of centrality c for pt — 4 GeV/c in black lines 
and pt = 6 GeV/c in red line. Dash-dotted line is for p a — 2 
GeV/c and dashed line for p a — 4 GeV/c. The solid lines 
are for the yields integrated over p a from 2 to 4 GeV/c. The 
data are from Ref. gtj for 3 < p^ is < 4 GeV/c and pf soc > 2 
GeV/c. 

To see the dependence on p a for fixed pt we show the 
yield in Fig. 25 for two representative values of c at 0.05 
and 0.35. These distributions are very similar to YJ^ 31 
in Fig. 3, which is the yield for fixed j3L = 2.9. Thus the 
result is the same whether we fix c or [3L. The inverse 
slope T a is therefore essentially what is shown in Fig. 4 
already. In Fig. 25 we have included two data points 
from Ref. (4lj . where recent results on near-side correla- 
tions have been reported. The data for central Au+Au 
collisions (0-10%) at 200 GeV are given for integrated jet 
yield per trigger (for — 1 < A</> < 1) with 3 < p^ lg < 6 
GeV/c and 1.5 < pf? soc < p^ E . Since our calculation is 
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FIG. 25: Yield per trigger in the near-side jet as functions of 
p a for two values of c and three values of p t . Data points are 
from Ref. [4l[; see text for details. 



for per-trigger yield of particles in the near-side jet av- 
eraged over all A(j> in the jet, we have divided the data 
by 2 (the range of A<fi), and include only the points at 
passoc > 2 GeV/c in the upper panel of Fig. 25. Our curve 
for pt = 4 GeV/c agrees very well with those two data 
points, which are averaged over the range 3 < pt < 6 
GeV/c. 
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FIG. 26: (Color online) A 3D plot of Y™* r (p t , p a , c) for c = 
0.05 (red), c = 0.35 (yellow) lowered by a factor of 10 , and 
c = 0.86 (blue) lowered by 10~ 2 . 



increase in the magnitude is seen when c is raised from 
0.05 to 0.35. Thus when the nuclear overlap is smaller, 
it is easier for the recoil jet to reach the away side and 
to produce a particle at pb- The shape of the pb dis- 
tribution is basically independent of centrality, since the 
hadronization process does not change with c. Figure 28 
shows a 3D plot of Y™ ay (p t ,pb,c), again with c = 0.35 
(in yellow) and c = 0.86 (in blue) lowered by factors of 
10" 1 and 10~ 2 , respectively. The near independence on 
p t is evident, while the increase with c is only from 0.05 
to 0.35, but not from 0.35 to 0.86. 
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FIG. 27: Yield per trigger in the away-side jet plotted in the 
same format as in Fig. 25 with p a replaced by p b . 




An overall view of Y£° ar (p t , p a , c) as a function of both 
pt and p a for three illustrative values of c is shown in 
Fig. 26. For clarity's sake we have multiplied the yield 
for c = 0.35 (in yellow) by 10 _1 and for c = 0.86 (in blue) 
by 10 -2 . The increase with p t is perceptible, while the 
dependence on c is negligible. 

For the away-side yield we use Eq. (|3"8"|) and obtain 
the results shown in Fig. 27, where a factor of about 2 



FIG. 28: (Color online) Same as in Fig. 26 but for 

Y:r y (Pt, Pb ,c). 

For fixed c and varying combinations of pt and pb, we 
can determine a quantitative measure of the antitrigger 
bias by calculating the average ((3t'}, where t' denotes the 
distance from the hard-scattering point to the away-side 
surface. In the calculation we identify j3t' with ln(k/q') 
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by virtue of G(q',k,L-t) inEq. ([22]). Figure 29 shows the 
results for c = 0.05 and 0.35. For fixed p t , (fit'} decreases 
with increasing pb as the scattering point is pulled closer 
to the away side. For fixed pb, that point moves closer 
to the near side, as pt increases, thus increasing {(3t'). 
The whole set of curves are lower at higher c. Thus Fig. 
29 provides a good description of the antitrigger bias. 
Note that the magnitude of (fit 1 ) is not large, less than 
0.5 even for c = 0.05. It is much smaller than the value 
f3L = 2.9 in Eq. (f2"5| . which is for the single-particle 
inclusive distribution. Again, we see that when a particle 
on the away side is required, the scattering point cannot 
be too far from the surface of the away side. At p t — 
Pb = 4 GeV/c, we have (f3f) « 0.2 for both values of 
c. That is just the value of {(3t) at pr — 4 GeV/c in 
Fig. 7, consistent with the scattering point being midway 
between the two sides. As we have learned from the 
preceding section, when p t = Pb the two jets produced 
are tangential jets near the rim of the overlap. As p t is 
increased, the scattering point can be embedded deeper 
in the interior, so {(3f) increases, but not very much. 
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FIG. 29: Average value of fit' where t' is the distance between 
the hard-scattering point and the away-side surface for various 
values of c,pt and pb- 

We can also present a more explicit description of the 
antitrigger bias in terms of parton momenta. We show in 
Fig. 30 the average momentum (k'}(p t ,pb,c) of the hard 
parton directed away from the trigger as a function of 
Pb for fixed c and pt- The increase with pb is now much 
slower than the increases with p {= Pt = Pb) in Fig. 
21. The difference between c = 0.05 and 0.35 is minor, 
as in Fig. 21. Focusing oiip t — A GeV/c, we see that 
(k 1 ) increases from ~ 8 to ~ 13 GeV/c, as pb increases 
from 2 to 6 GeV/c, the magnitude being significantly 
lower than the corresponding (k 1 ) in Fig. 11 for fixed 
/3L = 2.9. However, compared to (k) on the near side 
in Fig. 6, where (k) 8 GeV/c at p T = 4 GeV/c, (k') 
starts from about the same value at low pb but increasing 
persistently with pb, although p t is fixed. That is, (k') 
is always greater than (k) at fixed c and pt for one of 



the following two reasons. If pb is less than p t , then the 
longer path length on the away side due to antitrigger 
bias leads to higher (A/) despite the momentum balance 
k! = k in every hard scattering event. If pb is more than 
Pt, then clearly the jet momentum on the away side must 
on average be higher than on the trigger side. 




FIG. 30: Average value of the initial hard parton momentum 
directed at the away side for various hadron momenta pt and 
Pb for two values of centrality. 



^ lb 
o 

> 12 

CD 



8 
4 

s ^12 
9 8 



Q. 



c=0.05 




- p =4 GeV/c 
p!=6 , 



c=0.35 




3 p b (Giv/c) 5 



FIG. 31: Average value of the parton momentum at the away 
side surface for various p t , pb, and c as in Fig. 30. 

The behavior of (q')(pt,Pb,c) at the away-side surface, 
shown in Fig. 31, differs from that in Fig. 22 in the 
same way that Fig. 30 differs from Fig. 21. Compared 
to (k')(p t ,pb,c), the magnitude of (q')(p t ,pb,c) are, of 
course, lower, but the pt and pb dependencies are simi- 
lar. More revealing is the ratio (q 1 ) / (k')(p t ,pb, c) in Fig. 
32, which shows the effect of energy loss that decreases 
(difference from 1) with increasing pb due to decreasing 
path length, but increases with increasing p t due to in- 
creasing path length. The push-and-pull effect of pt and 
Pb is now clearly depicted in Figs. 30-32 that could not 
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FIG. 32: The ratio {q')/{k') for p t , p 6 and c as in Figs. 30 and 
31. 

We have been able to exhibit these characteristics of 
the medium effects by calculating theoretical quantities, 
such as ({3t'), (k') and (q'), for various values of p t , pb 
and c that are experimentally measurable. Beside the 
per-trigger yields Yj^ aI (p t ,p a , c) and Y™ ay (p t ,p b , c) that 
can be directly checked by experiment, it is possible that 
(q) and (q'} can indirectly be estimated by studying the 
total transverse momenta of charged particles in the near- 
and away-side jets. 



VIII. CONCLUSION 

We have made an exhaustive investigation of the prop- 
erties of dihadron correlation in jets produced in heavy- 
ion collisions. Our treatment is based on a reliable de- 
scription of hadronization through recombination on one 
hand, and on a realistic accounting of the medium ef- 
fect on the other. Only two free parameters are used to 
specify the distribution of dynamical path length £ for 
any centrality, and they are determined by fitting over 
100 data points of the n° inclusive spectra. All other 
parameters have been fixed by previous work in the re- 
combination model. Thus our study of the hadronic cor- 
relation with high-px trigger has very little freedom for 
adjustment, and for that reason we have been able to 
calculate unambiguously many quantities that reveal the 
medium effect on partons and the relationship among 
the momenta of hadrons that they produce. Not all the 
quantities calculated can be measured, but among those 
that can be checked by experiments, each encounter with 
existing data exhibits good agreement. It is therefore 
reasonable to conclude that the theoretical framework 
presented here offers a reliable description of one- and 
two-jet production at RHIC energy. 

One outcome of this study is the determination of the 
probability P(£, c) of a hard parton having a dynami- 



cal path length £ in a collision at centrality c, where £ 
plays the role of the suppression parameter f3L, except 
that it varies among all possible trajectories and density- 
dependent energy-loss factors, while j3L is the average 
over all £. The behavior of P(£, c) exhibited in Fig. 19 
for various c may be regarded as the fruit of the program 
to learn about the medium effect from heavy-ion colli- 
sions. Our calculation of the dihadron correlation in jets 
for each centrality would not have been reliable enough 
to compare favorably with data as in Figs. 24 and 25, if 
we did not have P(£, c) to link theory with experiment. 

Our study has shed light on the trigger and antitrig- 
ger biases, which are brief terms, referring to the com- 
plex issues involving nuclear geometry and medium sup- 
pression, that we have made more precise by examining 
the average momenta (fc), (k'}, (q) and (q') under differ- 
ent conditions. The ratios of those quantities reveal the 
different suppression factors on partons propagating to- 
ward the near and away sides. When the trigger momen- 
tum and the away-side associated-particle momentum are 
equal, we have learned, somewhat by surprise, that the 
yield is dominated by tangential jets, essentially indepen- 
dent of centrality. Revealing results such as that await 
precise verification by experiments. 

We have restricted our attention in this paper to the 
dependence on pt only, and to the range of pt large 
enough to leave out the consideration of ridge on the 
near side and of double humps on the away side. Ridge- 
ology is a separate subject in its own right and is treated 
elsewhere 28]. Lowering the pt range would contami- 
nate the trigger with medium partons through thermal- 
thermal recombination, a situation well within the ca- 
pability of our formalism to handle, but not considered 
here. Our treatment can also be generalized to include 
the azimuthal angle </>. The dependence of jet produc- 
tion on the trigger angle <p s relative to the reaction plane 
would be very interesting to study, as has already been 
initiated in data analysis by STAR [l^]. The widths of 
the jet peaks on the two sides, as well as the possible 
misalignment of the back-to-back jets, are challenging 
problems still to be investigated. 

A significant qualitative conclusion that can be drawn 
from the many ways of posing our questions on the 
medium effect is that it is hard experimentally to probe 
the interior of the collision zone. The profile of the nu- 
clear overlap is like a small island whose inhabitants may 
live uniformly throughout the island, but only the ones 
near the shore can react easily to activities in the sea. 
The detected trigger jets are produced mainly along the 
near side, like the coastal inhabitants responding to a call 
from the sea. If an away-side jet is detected at the same 
time as the trigger jet, the production point is moved 
along the rim so that the back-to-back jets are domi- 
nantly tangential jets, just as in the analogy where the 
inhabitants that can respond to calls from both sides are 
the ones on the coast with unobstructed vision of the 
two sides. In that sense jet tomography fails to probe 
the interior of medium, since rim production overwhelms 
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any signal arising from the interior in one- and two-jet 
events. The situation naturally suggests that three-jet 
events may reach the interior, analogous to the inland 
dwellers of the island being able to have equal, though 
harder, access to all points at sea. In e + e~ annihila- 
tion gluon jets were discovered in events where three jets 
were produced, each being at nearly equal azimuthal an- 
gles from the other two. Similarly, three jets originating 
from a common point in heavy-ion collision are possi- 
ble, such as in the gluonic process g + g — ► g + g + g. 
They are, of course, less abundant than two jets for the 
dual reason of higher order and enhanced suppression. 
However, the more serious experimental difficulty is to 
distinguish the 3-jet events from the background con- 
sisting of the double-hump structure in the conventional 
away-side distribution. The latter is due to TT re- 
combination, if the double humps can be related to the 
Mach cone, whereas the 3-jet event structure is associ- 
ated with shower partons in each jet. The subject is rich 
and worthy of attention, if the medium interior is to be 
probed by jet tomography. 

Another area of extension from this work is obviously 
heavy quark physics. In our sum over all parton types i 
in Eqs. (fT8|) . (fT9j) and (|23|) . we have limited to g, u, d and 
s, and their antiquarks. To consider charm quark, for ex- 
ample, we must redo everything from the beginning, in- 
cluding the shower parton distribution and single-hadron 
(D) spectrum. The basic formalism is, however, the same 



as we have given here. A new P(£,c) would have to be 
found for heavy quarks, and different possibilities of dijct 
correlation may reveal different medium effects [42[ . 

Finally, it is worth commenting that dijet correlation 
will undoubtedly be drastically different at LHC where 
not only numerous jets with px < 20 GeV/c will be cre- 
ated, but also the recombination of shower partons from 
neighboring jets can totally change the correlation be- 
tween hadrons [43| . The formalism used in this work will 
have to be thoroughly revised. Even the notion of jets 
distinguishable from background will require reexamina- 
tion. The challenging work ahead will, however, not be 
daunting but stimulating, both theoretically and experi- 
mentally. 
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